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The boundaries of the chart of nuclides contain exotic isotopes that possess extreme proton-to-neutron 
asymmetries. Here we report on strong evidence of °N, one of the most exotic proton-rich isotopes where 
more than one half of its constitute nucleons are unbound. With seven protons and two neutrons, this 
extremely proton-rich system would represent the first-known example of a ground-state five-proton 
emitter. The invariant-mass spectrum of its decay products can be fit with two peaks whose energies are 
consistent with the theoretical predictions of an open-quantum-system approach; however, we cannot rule 
out the possibility that only a single resonancelike peak is present in the spectrum. 
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Nuclei with large imbalances between their constituent 
numbers of protons and neutrons can have exotic proper- 
ties. The largest imbalances occur beyond the proton and 
neutron drip lines where the nuclear ground states (g.s.) 
are unbound. Because of the odd-even staggering of the 
drip lines induced by the nucleonic pairing, the shedding 
of unbound protons is usually terminated in an even-Z, 
particle-bound residue. Thus, just beyond the proton 
drip line, one is likely to find single-proton emitters for 
odd-Z isotopes and two-proton (2p) emitters for even-Z 
isotopes [1,2]. Even further removed, one finds 3p and 4p 
emitters. Presently, ’B, '7F, Na, and 3!K are the known 3p 
emitters [3-6] and °C and !8Mg have been observed to 
decay by emission of four protons [7,8]. 

Moving outward past the drip lines, the decay widths of 
the low-lying states increase, eventually melting into an 
unresolvable continuum as their lifetimes become commen- 
surate with typical reaction and single-particle timescales. 
Here, the very notion of the nuclear state becomes ques- 
tionable as the timescales are too short to talk about the 
existence of a nucleus. Indeed, as discussed in Ref. [9], if a 
collection of nucleons survives for less than about 1072? s, it 
should not be considered a nucleus. In this regime, the decay 
properties manifest themselves as scattering features rather 
than well-defined resonances. The maximum decay width at 
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the boundary for A x 8, based on single-particle timescales, 
is of order T + 3.5 MeV [10]. 

The 4p emitter 8C has one half of its nucleons in the 
continuum and the remainder constitutes an a cluster. The 
nucleus °N is even more extreme with an additional proton in 
the continuum. In nuclei such as 8C and °N, we might 
anticipate difficulty in classifying nuclear properties in terms 
of localized nuclear states. On the other hand, new insights 
can be gained from such unbound clusters of nucleons. In 
this Letter, we make a careful distinction between reso- 
nances (sharp peaks in the experimental scattering cross 
section) and resonant states (complex-momentum poles of 
the scattering matrix). For the theoretical classification of 
resonant states, see Fig. 1. The decaying resonant states in 
the fourth quadrant of the complex-k plane, lying close to the 
real k axis and having a real energy Re(E) > 0 and a width 
T = —2Im(E) > 0, can be interpreted as narrow resonances 
seen in experiments. The poles with Re(Z) < 0 and T > 0 
can be associated with subthreshold resonant states 
(0 > 45°). The antibound states have Re(E) < 0 andT = 0. 

Unlike the g.s. of 8C, the dineutron is not a genuine 
resonance (energy E > 0, decay width T > 0), but an 
antibound state (or a scattering feature) [11,12] located on 
the second Riemann energy sheet. Here the attractive 
interaction between the two neutrons is just insufficient 
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FIG. 1. Classification of resonant states based on the location of 
the scattering-matrix poles in the complex-k plane. The same 
arbitrary units (a.u.) are used on both axes. Bound, antibound, 
decaying, and capturing resonant states are marked, as well as 
narrow resonances (nr), broad resonant states (br), and subthresh- 
old resonant states (sr). The distribution of poles is symmetric with 
respect to the imaginary k axis because of time reversal symmetry; 
thus, capturing states are presented as the time-reversed decaying 
states. The dashed —45° line separates decaying resonant states 
from subthreshold poles. The inset shows the polar angle @ of the 
resonant pole in the complex momentum plane. 


to produce a bound state, but its nearly bound nature is 
manifested by enhanced n+ scattering just above 
threshold and by strong final-state effects. The diproton 
is a subthreshold resonant state (Æ < 0,P > 0) [13,14] 
and located below the —45° line in the complex-k plane 
(see Fig. 1). It is formally neither an antibound state nor a 
resonance (though closer to the latter), and again mani- 
fests itself by enhanced scattering strength and final-state 
effects. While there have been some recent suggestions of 
a tetraneutron resonancelike structure, what was observed 
[15,16] may instead be a final-state effect [17,18]. Similar 
situations might occur in the g.s. of °N and °He as 
presented in this Letter. 

The nucleus °N has three neutrons less than the lightest 
particle-bound nitrogen isotope !?N and one more proton 
than °C into which it decays. The neighboring isotope '°N 
has only been observed in three studies [19-21]. It has low- 
lying states which are single-proton resonances although 
their structure is not well established. Some indication as to 
the structure of °N can be gleaned from its mirror partner 
°He for which low-energy resonancelike states decay by the 
n + *He channel. 

Particular interest in °He is due to the parity inversion of 
the ground-state spin for odd N = 7 isotones with large 
neutron excesses when the second s;/7 neutron single- 
particle orbital intrudes into the p shell [22]. Most studies 
agree that there is a 1/27 resonance around 1.2 MeV above 
the n + *He threshold (see Ref. [23]), but there is less 
agreement about its width [23-25]. A number of studies 


find some 1/27 strength below this resonance [23,25—28] 
although this strength in some studies is not sufficient to 
justify a notion of a state [25,28,29]. 

Experiment.—An E/A = 69.5 MeV secondary beam of 
BO (4x 10° pps, purity 80%) was produced from the 
Coupled Cyclotron Facility at the National Supercon- 
ducting Cyclotron Laboratory at Michigan State Univer- 
sity. Charged particles created in the interaction with a 
1-mm-thick °Be target were detected in the high resolution 
array (HiRA) [30] consisting of 14 E-— AE telescopes 
covering scattering angles from 2.1° to 12.4°. The location 
of the HiRA telescopes was chosen based on Monte Carlo 
simulations which optimized the detection of low-lying 
resonances produced in projectile-fragmentation reactions. 
More details of the experiment and the simulations can be 
found in the Supplemental Material [31]. States in °N were 
produced in fragmentation reactions where three neutrons 
and one proton in total are removed from the projectile and 
the resonance is identified using the invariant-mass tech- 
nique. Data from this experiment pertaining to the first 
observation of ''O [49,50] and 'F [4] as well as other 
previously known isotopes [21,51,52] have already been 
published. 

Invariant-mass spectra.—The decay-energy (Qs,,) dis- 
tribution for all detected 5p + a events [Fig. 2(a)] is quite 
wide and contains no prominent “narrow” peaks, but there 
can be contributions from the decay of very wide °N states 
and nonresonant breakup. These events are produced in 
projectile-fragmentation reactions as their center-of-mass 
velocities are close to the beam value [53]. Given that 
the lowest-energy states in the mirror system °He decay to 
the n + *He(g.s.) channel, it is thus interesting to look 
at the mirror p + 8C(g.s.) channel. For each 5p + a event, 
we can remove one of the five protons in turn to create the 
five possible 4p + a subevents. The distribution of êC 
decay energy (Q,,,) from these subevents [Fig. 2(b)] has a 
peak corresponding to the g.s. of 8C at the expected value 
of Q4,. For comparison, the dotted green curve shows the 
shape of the °C peak obtained from fitting the more 
numerous 4p + a events (see Supplemental Material [31]). 

Possible °N > p + °C(g.s.) decay events were selected 
by requiring that at least one of the 4p + a subevents lies in 
the gate around the ground-state peak [Fig. 2(b)]. Based on 
the presented fit to the ground-state peak in this figure, 20% 
of the 4p + a subevents in this gate are associated with the 
background component. However, the gate on p + *C(g.s.) 
is cleaner than this as some of the events contribute to both 
the peak and this background through their different 
subevents. Indeed, 11% of selected events have two 
subevents in this gate which accounts for 1/2 of the 
background contamination. These numerical partitions 
are confirmed by our simulations. 

The final distribution, with the 8C(g.s.) gate, is largely 
restricted to Qs, < 10 MeV [data points Fig. 2(a)] and 
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FIG. 2. Invariant-mass spectra for °N and 8C. (a) Histogram 
shows the distribution of the decay energy Qs, obtained from 
5p +a events with the invariant-mass method. The data points 
show the distribution where an °C(g.s.) intermediate state was 
identified. (b) The data points show the distribution of C decay 
energy (Q4p) from the five 4p + a subevents in each 5p + a 
event. The solid red curve shows a fit to this distribution using the 
experimental °C line shape obtained from detected 4p + a events 
(dotted green curve) plus a smooth background (dashed blue 
curve). The gate used to select p + 8C(g.s.) events is indicated. 


appears to have two peaks. Hence it may be a doublet rather 
than a singlet although the statistics are marginal for this 
distinction. 

Theoretical models.—Since continuum effects are strong 
for both C and °N, we used the complex-energy Gamow 
shell model (GSM) to predict the location of the nuclear 
states of interest as it has been used to study many weakly 
bound and unbound systems [54—56] including the g.s. of 
8C [32]. GSM differs from the traditional closed-quantum- 
system nuclear shell model as it allows for bound, scatter- 
ing, and Gamow resonant states with outgoing asymptotic 
behavior to be treated on equal footing by implementing the 
Berggren basis [57]. 

Calculations for °C and °N were performed by assuming 
an a core surrounded by four and five valence protons, 
respectively. We used the same valence-space Hamiltonian 
as described in previous GSM studies [32,33,58] with the 
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FIG. 3. Level diagrams of ŝC and °N obtained experimentally 


and calculated with the GSM. Energies are given relative to the 
4He threshold. The level diagram of °He, the mirror partner of oN, 
is shown in the inset. The 1/2+ antibound state in °He is shown 
with a wavy line to indicate its status more as a scattering feature 
rather than a real state. The proposed 1/2* state in °N is shown 
with both straight and wavy lines to indicate the uncertainty as to 
its nature (a resonance or scattering feature) while the GSM 
interprets it as a broad resonant state. 


parameters given in the Supplemental Material [31]. The 
maximum number of particles allowed in the continuum is 
2 for the results presented; it has been checked that more 
particles in the continuum barely change the conclusion for 
the A = 8 system. The predicted Q4, values for the g.s. and 
the 2+ first excited state of °C shown in Fig. 3 with the 
relevant levels for their decay are around 3.7 and 8.3 MeV, 
respectively (see Supplemental Material [31]; see also 
Ref. [59] for complex-scaling predictions of °C). 

The properties of the low-lying states in °N were 
cross-checked by the Gamow-coupled-channel (GCC) 
method [34,60], which is a complex framework utilizing 
the Berggren basis. Based on large proton-decaying 
spectroscopic factors Con = 0.87 and Shin = 0.80 
according to GSM calculations), °N can be described as 
a ŝC + p two-body system. 

Discussion—While we have cleanly selected 
p +*C(g.s.) events, the extra proton beyond the °C decay 
products can result from either °N resonant decay or 
alternatively is produced directly in the first step of the 
projectile-fragmentation reaction. To form °N, one proton 
and three neutrons need to be removed from the projectile 
and thus one prompt proton would be expected unless it is 
bound in a d or t cluster with some of the promptly removed 
neutrons. Also, there can be contributions from events 
where the projectile promptly breaks up into the channel 
2p + 3n + ®C(g.s.). While most of the promptly emitted 
protons are expected at larger angles than sampled with 
HiRA, we do expect some background from these prompt 
particles [53]. 
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The background from all these sources of prompt 
protons has been extensively studied in Ref. [53] for this 
and other resonant-decay channels associated with this 
projectile. The Qs, dependence of this background can 
be obtained from mixing p and °C decay products from 
different p + *C(g.s.) events appropriately weighted to 
include correlations between them. However, this predicted 
Qs, distribution is too wide to explain all the observed 
events and thus substantial contributions from °N reso- 
nances are required [53]. Indeed, the possibility that 
statistical fluctuations in this background can explain the 
observed spectrum can be ruled out with a significance of 
roughly 5ø, i.e., at the discovery threshold for a new 
resonance (see Supplemental Material [31]). 

The GSM predicts 1/2* and 1/27 states in the region 
below Qs, = 10 MeV where the p+ °C(g.s.) yield is 
strongest. The predicted ground state has the same parity 
inversion as observed for !'Be [22]. With two states and the 
magnitude of the background to vary, there are seven fit 
parameters, which is about the same as the number of data 
points in the peak region. Thus the experimental statistics 
are presently not sufficient to fully constrain the location of 
both resonance poles (Qı, — iI'/2) and determined their 
nature (genuine or subthreshold resonances) from Fig. 1. In 
the following, in order to show some possible inter- 
pretations of the data, we present some constrained fits 
which avoid overfitting. 

At one extreme, we can assume just a single peak plus 
background. Using R-matrix line shapes [35] for the p + 
8C(g.s.) channel, the qualities of the fits were very similar 
for both s and p resonances. The resonance pole for a single 
s-wave fit is at Qi, = 1.22(16) MeV [Q;, = 4.70(16)] 
and T = 2.59(23) MeV. Such a state is approaching the 
diffuse borderland between a broad resonance and a 
scattering feature. However, this solution is disfavored as 
the magnitude of the background has been reduced to zero 
to achieve the best fit. The result is not much better for a 
single p-wave resonance [Fig. 4(a)] with the resonance pole 
given by Q), = 2.01(16) MeV [Qs, =5.49(16) MeV] 
and I = 2.28(23) MeV. In this case, the fitted background 
(dashed blue curve) contributes only 14% of the total yield. 
As other resonant channels from this 'O-induced re- 
action had backgrounds of magnitude 30%-54% [53], these 
single-peak solutions are disfavored. 

Figure 4(b) shows an example of a possible two-peak 
fit assuming real p + %C(g.s.) resonances described by 
R-matrix line shapes. Here, the R-matrix parameters of 
the 1/27 peak are fixed so that its pole is consistent with the 
GSM prediction. With this constraint, the fitted values for 
the 1/2* strength are Qi, = 2.75(21) MeV [Q5, = 
6.23(21) MeV] and T = 0.58(44) MeV, which can be 
compared to the GSM predictions of Q5, = 5.56 MeV 
and I’ = 1.74 MeV. The fitted background from prompt 
protons (dashed blue curve) now explains the observed 
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FIG. 4. Fits (solid red curves) to the selected p + *C(g.s.) 
distribution. (a) Fit with a single peak parametrized by an ¢ = 1 
R-matrix line shape. (b),(c) Two peak fits where the contributions 
from both levels are given by the dotted green curves. The fit in 
(b) was using a p + °C(g.s.) R-matrix line shape where energy 
and width of the upper peak are set to the GSM predictions for the 
J” = 1/27 state of °N. The dependence of the detection effi- 
ciency as a function of Qs, as determined in our Monte Carlo 
simulations (see Supplemental Material [31]) is shown by the 
dashed magenta curve. The fit in (c) was obtained where the fixed 
width of the 1/27 was halved from its value in (b) and the 1/2* 
strength is described by a subthreshold resonant state using GCC 
line shape. In all panels, the p + 8C(g.s.) threshold is indicated 
by an arrow and the fitted background from nonresonant protons 
is shown by the dashed blue curves. 


yields below Q5, = 5 MeV and the high-energy tail above 
10 MeV. This background accounts for ~40% of the yield, 
which is similar to values of 30%-54% obtained in fitting 
the 11N = p+°!°C invariant-mass distributions [53]. 
While this fitted width of the 1/2* state is smaller than 
the GSM prediction, it can be increased by decreasing the 
fixed width of the 1/27 state. For example, Fig. 4(c) shows a 
fit where the intrinsic width of the 1/27 state has been 
halved, but now the 1/2* strength is described by the GCC 
line shape as a subthreshold resonant state. 

To provide insights into the nature of the 1/2* state in °N, 
GSM calculations were performed for its mirror partner, 
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°He, using the same Hamiltonian parameters. Figure 3 
shows that the experimental spectrum is reasonably repro- 
duced, although the prediction for the 1/27 state is some- 
what too high in energy. The previous studies using the 
Berggren basis have indicated that the 1/2* state in °He is a 
resonance [33,36], in which some bound poles were 
introduced in the basis to stabilize their results (see 
Supplemental Material for details [31]). However, when 
the continuum effect is properly taken into account in the 
current work with a deformed scattering contour [37], one 
can generate an antibound 1/2+ pole in °He, an approach 
yielding the best agreement with experimental data. 

Although it is not possible for antibound poles to 
emerge in proton-rich nuclei due to the Coulomb inter- 
action [34], we followed the same procedure as in °He to 
determine if the 1/2* state in °N is a broad resonant state 
or subthreshold resonant state. The predicted 1/27 state 
has an energy of E = 2.08 MeV above the 8C + p thresh- 
old and I = 1.74 MeV, which indicates a broad resonant 
state as defined in Fig. 1. It must be noted, however, that the 
1/2* state in GSM is very fragile with respect to changes of 
the Berggren basis due to the truncations (particle- 
hole, discretization of the continuum) employed (see 
Supplemental Material for details [31]), so we cannot with 
certainty rule out a subthreshold resonant state. 

Both °C and °N shed their excess protons by sequential 
steps of either single-proton or prompt 2p decays. The 
4p-emitter '*Mg also decays in this manner [8], suggesting 
this behavior is typical for isotopes at the limits of 
existence. Our analysis suggests that some observed 
structures, sometimes interpreted in terms of nuclear states, 
should be viewed rather as fleeting features lying outside 
the chart of nuclides. 

Conclusions.—We have found strong evidence for the 
exotic nuclide °N produced in the fragmentation of a 13O 
beam. The invariant-mass spectrum of detected 5p + a 
events, each containing an *C(g.s.) intermediate state, 
contains a structure which cannot be explained at the 
So level by statistical fluctuations of the expected back- 
ground from other coincident protons liberated in the 
fragmentation event. This observed structure can be inter- 
preted as two °N peaks, although a single-peak solution 
cannot be totally discounted. The °N state(s) would be first- 
known case(s) of 5p emission from ground-state and 
low-lying resonances. 

This nuclide has also been studied theoretically in 
the Gamow shell model where the important effects of 
the continuum are included. The predicted locations of the 
1/2* and 1/27 resonant states are in excellent agreement 
with the location of the observed structure giving further 
evidence for the preferred two-peak solution. In the model, 
the 1/2+ ground state of °N is the mirror of an antibound 
state in °He, and is most likely a broad resonant state rather 


than a subthreshold resonant state. However, the latter 
cannot be completely ruled out in both the experiment and 
the theory. 
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